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a^ ■ Abstract. The long-time behaviors of the velocity autocorrelation function (VACF) 

! for sheared fluids are investigated theoretically and numerically. It is found the 

existence of the cross-overs of VACF from t^'^^^ to i^'' in sheared fluids of elastic 
particles without any thermostat, and from i^''/^ to i^('*+2)/2 jn both sheared fluids 
of elastic particles with a thermostat and sheared granular fluids, where d is the 
C^ ' spatial dimension. The validity of the predictions has been confirmed by our numerical 

ff , simulations. 
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The slow relaxation of the current autocorrelation functions to the equilibrium 
state is one of the most important characteristics in nonequilirbium statistical physics 
[lj-[5j. It is known that the long-time tails are caused by the correlated collisions 
which cause anomalous transport behaviors such as the system size dependence of the 
transport coefficients P, [7j. The analysis of such behaviors is still an active subject in 
nonequilibrium physics. Although there is the long-time tail t~'^^'^ with the time t and 
the spatial dimension d even at equilibrium, the tail under an external force such as 
a steady shear has different feature from that at equilibrium. Indeed, we believe that 
the viscosity rj depends on the absolute value of the shear rate 7 as 77 — ?7o = rj'log^j) 
for d = 2 and i] — 'Qo = r]"'j^^'^ for rf = 3 in the sheared ordinary fluids at a constant 
temperature, where tjq is the Newtonian viscosity at equilibrium, and rj' and rj" are 
proportional constants [HI [9]. These results are obtained from the assumption that the 
current autocorrelation functions for t > 7"^ decay much faster than the conventional 
tail. Nevertheless, because of the lack of systematic studies in the long time region, we 
arise the following question: (i) What are the long-time behaviors of the autocorrelation 
functions for t > 7""^? 

The system becomes unsteady because of the increment of the temperature due 
to the viscous heating effect when we add the shear to a system consisting of elastic 
particles. We, then, sometimes introduce a thermostat to keep an isothermal steady 
state of the fluid under the shear [Tn]-[l3]. We, however, do not know the details of 
roles of thermostats in the sheared fluids. On the other hand, the sheared granular fluid 
in which there are inelastic collisions between particles can be regarded as one type of 
isothermal fluids, because the sheared granular fluid can keep a constant temperature 
under the balance between the shear and the dissipation due to inelastic collisions. 
Indeed, our recent paper has conflrmed that an unifled description of both sheared 
granular fluids and sheared isothermal fluids with a thermostat can be used for equal- 
time long-range correlations as long as the systems keep uniform shear flows [H]. We, 
thus, arise the second question: (ii) What are actual relations among the sheared fluid of 
elastic particles without any thermostats, the sheared isothermal fluid of elastic particles 
and the sheared granular fluid! 

The interest in the role of the long-tails in granular fluids is rapidly growing \T5\- 
[22] . In particular, Kumaran suggested the existence of an interesting fast decay of the 
autocorrelation functions, t~3d/2 ^^ sheared dense granular fluids [l5l[T6], while Otsuki 
and Hayakawa predicted that the correlation of the shear stress decays as t"('^+^)/2 [ T7] . 
These predictions are nontrivial, but the exponents for the tails are not conflrmed in 
the recent experiments [H] . On the other hand, the existence of conventional long-time 
tails in velocity autocorrelation function (VACF) has been confirmed, while the absence 
of long-tail in the correlation of the heat flux is found p-9j. Thus, we encounter the third 
question: (iii) What is the true long-time tail in the sheared granular fluids'! 

In this letter, to answer the above three questions we study the long-time behaviors 
of VACFs in sheared elastic particles with or without thermostats and sheared granular 
particles. Our theoretical method is based on the classical one developed by Ernst et 
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al. [21 lin] , and the theoretical results will be verified from our numerical simulations. 

Let us consider the system consists of A^ identical smooth and hard spherical 
particles with the mass m and the diameter a in the volume V. The position and 
the velocity of the i-th particle at time t are rj(t) and Vj(t), respectively. The particles 
collide instantaneously with each other with a restitution constant e which is equal to 
unity for elastic particles and less than unity for granular particles. When the particle i 
with the velocity Vj collides with the particle j with Vj, the post-collisional velocities v* 
and V* are respectively given by v* = Vj — |(l + e)(e-Vjj)e and v* = Vj + |(l + e)(e-Vjj)e, 
where e is the unit vector parallel to the relative position of the two colliding particles 
at contact, and Vjj = Vj — Vj. 

We assume that the velocity profile is given by Ca{r) = 'yy6a^x in our system, where 
the Greek suffix a denotes the Cartesian component. In this letter, we discuss the 
following situations: (a) A sheared system of elastic particles without any thermostat, 
(b) a sheared system of elastic particles with the velocity rescaling thermostat, and (c) 
a sheared granular system with the restitution constant e < 1. We abbreviate them 
the sheared heating (SH), the sheared thermostat (ST), and the sheared granular (SG) 
systems for later discussion. To avoid the consideration of the contribution from the 
potential terms, we assume that the sheared fiuid consists of a dilute gas of particles. 

We are interested in VACFs for d-dimensional systems [231 121! 

. N . d 

cl^\t)^-Y.Wr,MvUt)). c^'\t) = -Y^c^{t), (1) 

where we have introduced v'i^{t) = fj,a(0) — CQ,(rj(0)), and C^'^^t) is the superposition 
of C^\t). VACFs are expected to be represented by the fiuctuation of hydrodynamic 
fields around uniform shear fiow characterized by the uniform shear velocity CQ,(r) and 
the homogeneous temperature Th [21 [T9| as 

C^\t) ^ j dVMv'.K^^ j ^^uUt)P^^{tl (2) 

where tik.al^), Pv.it)^ Vq and /o(fo) are the Fourier transforms of a-component of the 
velocity field and the probability distribution function of a tracer particle, the initial 
peculiar velocity of a particle defined by Vq = vq — c(ro) with the initial velocity vq, 
and the initial velocity distribution function, respectively. Here, we should note that 
the wave number k is taken to satisfy the Lees-Edwards boundary condition [25]. The 
validity of eq. ([2]) in the nonequilibrium situations has been verified in Ref . [19] . Thus, 
we need to estimate the behaviors of Uk(t) and P]i{t) to obtain VACFs. 

As was shown in the previous studies [21 [Hj, we adopt the linearized hydrodynamics 
of fluctuating fields around the uniform shear fiow to evaluate VACFs. As in the 
equilibrium case, we decompose Uk(t) into the longitudinal mode Mk||(^) and the 
transverse mode u^\^{t) as Uk(t) = U]^\\{t)e\\ +I]f~^ u^'^{t)eY , where we have introduced 
eii = k(-7t)/fc(-7t), and ej^ = {e„, - (e,, • e^Je,, - Er'(e? ■ e„Je?}/Ar with 
k(7t) = k + ^tKey [25] and AT = |e„^ - (ey • e^Je,, - Ef^l^? • e^Je?! . Here, 
e^. is the unit vector parallel to aj direction. In this letter, we take ai as y. We 
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introduce w^j_{t) = u^\_{t) / J Tnit) , and the time evolution equation of Wkj_(t) is given 
by 



(*) i+\ — ,,* Fr(+\ui ^,+^2„„W 



a,<i(t) = - v*JTH{t)k{-W<iit)+inK-it).t), (3) 



where we have introduced F(k, t) = {g^ ■e.x)'U]<iy{t) / JTH{t)+5i^ikxk^/ ku\i\\{t) / JTH{t) + 
u*{ik:,{ef ■ey)+iky{ef ■e^)}Tk{t)/TH{t) with k± = Jk^ - k^ Here, T^{t) is the Fourier 



y 
transform of the temperature, and z/* is unimportant part of the kinetic viscosity which 

is independent of the temperature. 

We note that the transverse mode w^_^{t) is not independent of the fongitudinal 

modes for sheared fluids because of the existence of F(k, t). Nevertheless, we can ignore 

the mixing effect by assuming the small shear rate satisfying 7 <^ v* JTH{t)k{~'^t)'^ . It 

should be noted that the condition 7 <^ u* JTH{t)k{—'yt)'^ is not satisfied in the limit 



fc ^ 0. However, k has the infrared cutoff 2-k/L with the system size L for the sheared 
system. Thus, the condition can be used when we assume small enough 7 or high enough 
initial temperature Th{0) for SH and ST, and 1 — e^ <^ 1 for SG, as in Ref. p3]. 

We expect that the longitudinal mode proportional to ey is irrelevant for sheared 
fluids because of the existence of the sound wave in sheared fluids even when we consider 
SG. This is contrast to the case of freely cooling granular fluids [T9] . 

The time evolution equation of Pk{t) in sheared fluids is 

dtP^it) = - D*^/T^)k{-^tfP^{t), (4) 

where D* is unimportant part of the diffusion coefficient, which is independent of the 
temperature. 

From eqs. ([3]) and (jl]), it is obvious that the time evolution of the homogeneous 
temperature Tff(t) plays a key role in the long time behavior of VACF. When we 
consider a system of SH, the temperature increases by the viscous heating. In this 



case, the temperature obeys dTn/dt = 2mi'* ^/^H'^'^ / d and its solution is JTuit) = 
Th{0) + mu*'^'^t/d. On the other hand, Tnit) should be a constant for either ST or 



SG from the balance between the viscous heating and the energy dissipation. Therefore, 
we reach an important and an unexpected conjecture that the behavior of VACFs for 
nearly elastic SG is the same as that for ST, but is different from that for SH. The 
difference between nearly elastic SG and ST appears through the fact that the granular 
temperature in SG is determined by the shear rate as Th oc 7^, while the temperature 
in ST is basically independent of the shear rate. This conjecture is used in the analysis 
of the spatial correlations in SG |14j . 

Thus, the solution of eq. ([3]) with dropping F(k, t) is given by 

^kl(^) = wi^iiO)e-''V^^W^^^), (5) 

where 5(k,t) = Ai{t)k^ - A2{t)k^ky + A^{t)kl. Here, Ai{t) = t(l + (3t/2), A2(t) = 
7^2(1 + 2(3t/3) and A^it) = 'yHy3{l + 3f3t/4) with f3 = mz/*7V(d^Tf^(0)) for SH, and 
Ai{t) = t, A2{t) = 7^2 and Asit) = ^H^/3 for S T and SG. The solution of eq. (gD is 
similarly obtained as Pk(t) = Pk{7t)(0) exp[-L)* Jrff(0)5(k, t)]. 
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Substituting these solutions into eq. ([2]) with the aid of T^ = (m/c?) /(ivo(vo 
c)Vo(^o), Uk(0) ~ {V/N)vq and Pk(0) ~ 1, we obtain 
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^ J dkM^{k,t)e-^(^''\ (6) 



\ TniO) 



where we have introduced 

for a ^ y. Here, we note that the proportional constant in eq. ([6]) depends on the 
restitution coefficient e through z/*, D* and Tnit) in eqs. ([3]) and (jlj). Equations ([6]) 
and ([7j) describe VACFs for the sheared fluids in all-time region. 

In the short-time regime t < 7~^ VACFs, eq. (Q, are reduced to the known result 
for the fluid at equilibrium (EQ) C^\t) ~ t"*^/^ for all situations. On the other hand, 
the long time behaviors of C^\t) obtained from eq. ([6]) for t > 7"^, which depend on 
the situations, differ from those at EQ. 

Let us flrst consider the asymptotic behaviors of SH for t ^ 7^^, where B{t) in 
eq. IQ is dominated by A3{t)k'^ ~ j'^pt'^kl/A. To keep the contribution of this term 
we introduce the scaled wave number k' as k'^ = kx^/P'jt'^ and k'^ = ka^/Pt for a ^ x. 
Thus, VACFs for t ^ 7^^ approximately satisfy 

cf\t) ~ rH-\ cf\t) ~ r't-'. c('\t) ~ r't-V2 (s) 

for d = 2, while VACFs for d > 3 are approximately given by 

Next, let us consider the asymptotic behaviors of ST and nearly elastic SG for 
t ^ 7~^, where B{t) of eq. (Ej) is dominated by A3{t)kl ~ 'jt^kl/3. Similar to the case 
of SH, we introduce the scaled wave number k" as k^" = k^'jt^^'^ and kj' = k^t^f'^ for 
a ^ X. Thus, two-dimensional VACFs of ST and SG for t ^ 7"^ are identical to eq. 
(IH]), while VACFs for d > 3 are approximately given by 

C(f>5)^t) ^ C^''^'\t) ~ ^-lt-(rf+2)/2. (10) 

It should be noted that the long-time behaviors of C^f^{t) for SH, ST, and SG differ 
from those of C^\t) with a ^ y only in the case of rf = 2. To explain this different 
scaling of C^'^ (t) , we note that the scaled wave number k'^ or kx" differs from that of k'^ 
or ka" with a ^ x, and Mj^(k) given in eq. ([7]) is proportional to kj_ = k"^ — k^, which 
reduces to kj_ = k^ ior d = 2, and k"^ = k^ + J2a^x,y ^a for rf > 3. To demonstrate the 
difference, let us restrict our interest to ST. In the vicinity of k^ = 0, My(k) behaves 
as My{k) ~ A;i ~ {j:a^x,yka}/t for rf > 3, while My{k) ^ kl ^ ^^"7(7^^) for ^ = 2, 
which lead to the asymptotic behavior of Cjf-^^ ~ ^-it-{'i+2)/2 g^j^j C^\t) ~ 7~^t~^. 
Thus, we obtain the characteristic long-time behavior of C^^-* (t) among C^'^-* (t) . 

To verify the validity of our theoretical predictions we perform the simulations 
of hard spherical particle for d = 2 and 3. In our simulation, we flrst prepare the 
equilibrium state with the temperature T/ at t = to < 0. For t > to, we impose the 
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shear flow by the Lees- Edwards boundary condition to avoid the shear-band instabihty, 
and measure the velocity correlation from t = 0, where we have confirmed that the 
homogeneous shear flow with the temperature T^ is realized, m, a, and the temperature 



Tg at t = are set to be unity. Thus, the unit time scale is measured by Tg = aJm/Ts. 
In our simulation, the number of the particles we use is 65536 for the calculation of 
VACFs with the ensembles of 1600 different initial conditions for d = 2. The number of 
the particles we use is 262144 for the calculation of VACFs with the ensembles of 170 
different initial conditions for d = 3. We have already checked that the system size is 
large enough that any finite size effects are not observed. We adopt the area fraction 
z/ = z/c/2 with the closest packing fraction of particle Uc for d = 2 and u = z/c/1.9 for 
d = 3. The densities i^ we chose might be high, but our method can be extended to 
high density cases [Ml 
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Figure 1. (a) \ci^Ht)\ as the function of t for ST and EQ with d = 2. 
as the function of t for ST and EQ with d = 2. (c) C'^^(t) as the function of t for SH 
with d = 2. (d) C(2)(i) as the function of t for ST and SG with d = 2. (e) C'^^^t) 
as the function of t for ST and EQ with d = 3. Here, we use N = 65636, v — i^c/2 
and e = 0.99 for SG in (a)-(d), and N = 262144, i^ = 1/^/1.9 in (e). The shear rate is 
chosen as 7 = 0.2t-^ for SH, ST, and SG in (a)-(d), and 7 = 0.6t-^ for ST in (e). 



We use the parameters 7 = 0.2r~^ for d = 2, and 7 = 0.6t^^ for d = 3 in SH, ST, 
and nearly elastic SG. In the case of EQ, we use 7 = 0.0. For SG, we use the restitution 
coefficient e = 0.99 to set the temperature unity. Here we adopt the velocity rescaling 
method for ST where the velocity is rescaled for every time duration At/r^ = 0.01 to 
keep the constant temperature. 

Figures [1] (a) and (b) exhibit the numerical results of C^^-*(t) and C^^\t) for EQ and 
ST with d = 2. The theoretical curves in Figs. [11(a) and (b) for ST are obtained from eq. 
([6]) with two fitting parameters for the amplitude and the time scale. C!^^ (t) and C^^^ (t) 
for ST deviate obviously from those for EQ in the long time regime. These behaviors 
for ST are almost on the theoretical curves, which strongly supports the validity of our 
theory. We also stress that Fig. [1] (b) evidently supports the existence of the theoretical 
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cross-over of C^\t) from t~^ to t~^, though the numerical cross-over of C^\t) from t~^ 
to t~^ in Fig{T] (a) is obscure. From Fig. [1] (a), we confirm that |C^^''(t)| for ST rapidly 
decreases near t/Ts ~ 10 to become negative, while C^\t) is always positive, which is 
consistent with the theoretical prediction. 

Figures [^c)-(e) show the numerical results of C'''^-'(t). All the theoretical curves are 
obtained from eq. ([6]) with one fitting parameter for the amplitude. From Figs. [1] (c) 
and (d), we confirm that VACFs for SH, ST and SG with d = 2 are also consistent with 
the theoretical result. From Fig. [1] (e), we find that VACF for ST with d = 3 deviates 
from the line t~^/^ and VACF for EQ in the relatively long-time region. The data are on 
the theoretical curves again, which strongly support the validity of the theory, although 
the asymptotic behaviors for t ^ 7"^ could not be confirmed. 



One may be skeptical about the condition 7 ^ u* JTH{t)k{—'yt) to ignore 
F(k,t) in eq. ([3]), although the agreement between the theory and the simulation is 
good. For more precise analysis, one can solve the eigenvalue problem of linearized 
hydrodynamics around the uniform shear flow [U [T71 ^7\ . From the eigenvalue analysis, 
the transverse mode wl^\_{t) can be represented by the linear combination of two 
eigenvalues A"^(k(— 7t)) and A'^(k(— 7^)) and associated eigenvectors, where we have 
introduced A^(k) = iy*^/T^e - ^Kky/k'^ and A^(k) = v*^T^k'^. It should 

be noted that the transverse mode w\^l_{t) = W7^_j_(0) exp[— /q (isA"^(k(— 7s))] always 
exists, and w^!j (t) associated with A'^(k) does not exist for two dimensional case 
while w|;,*j (t) is d — 2 degenerated for d > 3. Since A'^(k) and its associated 
eigenvectors are identical to those discussed in eq. ([3]) with neglecting F(k,t), we 
should estimate the contribution from A"^(k) and w\^l{t). Substituting Pk(^) and 
w\^j_{t) with u\^l{t) = \/TH{t)wl^l{t) into eq. ([2]), we obtain its contribution to 



CW(t) as C^{t) oc Jdk{{k''-'itk,ky)/k{-^tf}exp[-J^ds^TH{s){u* + D*)k{-'is)% 
where we have used jQds'ykxky{'js)/k{'ysY = hi[k{'yt)/k]. C^it) obviously satisfies 
C^{t) ~ t^'^/^ for t <^ 7~^, while the long time asymptotic form for t ^ 7^^ is given 
by C^it) ~ j~H'(d-+'^)/^ for ST, where the proportional constant is J dK[{Ky + Kl - 



K.,Ky)/{{Ky - K,y + Kl}] exp[-^TH(0)(z/* + D*){K^ - K^Ky - 2^2/3)] for three 
dimensional case with K^. = k^i^jt^^"^) and Ka = kat^^'^ for a ^ x. Similarly, the 
asymptotic forms of C^{t) for SH and SG are proportional to those analyzed in the 
text. Therefore, we believe that the contributions from F(k, t) is only the change of the 
amplitude. The details of eigenvalue analysis will be discussed elsewhere [27] . 

Our result is contradicted with Kumaran's prediction [I5l |16]. This discrepancy 
does not disprove Kumaran's scaling, because ours is for dilute nearly elastic SG but his 
prediction is only valid for sheared dense granular fluids. However, we believe that our 
result is still valid even for dense sheared systems. Indeed, the universal results between 
SG and ST has been confirmed for the equal-time long-range correlation functions [9l[T4] 
even for dense systems, when the systems satisfy the uniform shear. We will discuss the 
unified description on both the long-range correlation and the long-tail elsewhere [27] . 

In conclusion, we have analytically calculated the behaviors of VACF in the sheared 
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fluids, which is represented by eq. (E])- (i) We have predicted the existence of the cross- 
over of VACF from t~'^/^ to t"*^ for sheared heating elastic particles (SH), while the 
cross-over is from t~'^/^ to t~('^+2)/2 for the sheared elastic particles with a thermostat 
(ST) and the nearly elastic sheared granular particles (SG) through the hydrodynamic 
analysis, (ii) We have also theoretically predicted that the behavior of VACF for ST is 
almost equivalent to that for SG, while SH is different from that of ST. (iii) As stated 
in (ii), we have confirmed that the hydrodynamic properties of dilute and nearly elastic 
SG are not unusual, which is contrast to Kumaran's prediction C{t) ~ t~^'^^'^. 
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